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It has been shown recently that negatively charged particle multiplic- 
ity distributions in e^e~ annihilation into hadrons [1-5] and in the lepton- 
nucleon scattering [6] are well described by the modified negative binomial 
distribution (MNBD). It has been shown also [7] that the MNBD and its 
simple extension EMNBD quite well describe charged particle multiplicity 
distributions in restricted (pseudo)rapidity intervals in e^e~ annihilation into 
hadrons and in e^p collisions at HERA energies. The aim of this paper is 
to show that the EMNBD describes negatively charged particle multiplicity 
distributions in e+e^ annihilation, measured in the whole phase space, as 
well as the MNBD. 

Let us remind that the MNBD can be defined by the probability gener- 
ating function 

where P„ is the probability to produce n particles, 

ip{x) = l-p{l-x) (2) 

and A;, A, p and r are parameters connected to the mean multiplicity < n > 
by the relation 

< n >= k p {r - A). (3) 

It has been assumed in the toy model proposed in [1,3,4] that the parameter 
k is the number of sources of particle production at some initial stage of 
the interaction; these sources develop independently of each other according 
to some branching process (characterized by the parameters r and A) and 
during the branching process intermediate neutral clusters are produced. 
The parameter p is equal to the cluster decay probability into a charged 
hadron pair and {1—p) is the probability of cluster decay into pair of neutral 
hadrons. The energy independence of the parameter A (or the product Ap), 
observed in the papers [1,3,4,5] can indicate that the branching process is 
the pure birth branching process and A can be fixed at —1. The probability 
generating function for the EMNBD has also the form (1) with ^{x) replaced 
by 

^2{x) = l-ei{l-x)-e2{l-x^), (4) 

where the parameters Si and 62 can be considered as the cluster decay proba- 
bilities into one and two pairs of charged hadrons respectively. The EMNBD 



transforms into the MNBD when £2 = 0, in this case ei = p. The probabih- 
ties Pn for the EMNBD can be calculated using the formulae given in [7]. 

The results of the EMNBD fits to the negatively charged particle multi- 
plicity distributions in e~^e~ annihilaion into hadrons [8-19] are given in the 
table 1. The parameter A was fixed at the value A = — 1 and the parameter 
r was calculated from the mean charged multiplicity < n > using the relation 

< n > , . 

r = A + — . (5) 

k (£1 + 252) ^ ^ 

The integer parameter k has been tested in the interval from 1 to 9, and the 
parameters ei and £2 have been assumed to be nonnegative. The errors for 
the parameter k were calculated using the quadratic interpolation for the x^ 
dependence on k on both sides from the Xmin- ^^^ can see from the table 1 
that the quality of the EMNBD fits is qood. The x^ for the EMNBD fits 
are in general smaller than the x^ for the MNBD fits [4,5], shown in the last 
column of the table 1. It is necessary to note that the x^ I^DF values for the 
fits should be considered just indicative, since the full covariance matrix is not 
given in the experimental publications and therefore the proper treatment of 
the correlations between measurements of the neighbour multiplicities is not 
possible. 

The energy dependence of the parameters Si and £2 is presented in the 
fig. 1. These parameters appear to be energy independent, if one excludes 
first and last energy points with ^i equal 3 and ~ 133 GeV. The big values 
62 at these energies can be explained by statistical fluctuations. One should 
note also that for these energies practically the same x^ jNDF values are 
obtained for the EMNBD fits with the parameters e\ and £2 fixed at the 
average values ~ 0.65 and ~ 0.16 respectively (not shown). 

The energy dependence of the parameter k for the EMNBD fits is com- 
pared in the fig. 2 with the energy dependence of the parameter k for the 
MNBD fits. The k for both parametrizations rise almost linearly with log(y^) 
at energies below ~ 30 GeV and seem to approach some asymptotic value 
~ 7 or 6 at higher energies. 

At the present level of the experimental precisions in the multiplicity 
measurements both the MNBD and EMNBD parametrizations look more 
or less the same, this is explained by the smallness of the £2 with respect 
to the £1. The good quality of fits is expected also for the next iteration 
when one adds the term £3(1 — x^\ responsible for the cluster decay into 



three charged hadron pairs, to the function ^p{x); indeed the probabihty £3 
is expected to be smaller than the 62- These iterations remind the Pade 
approximants (ratios of the polynomials) of the increasing order, used in the 
calculational mathematics for the function approximation. The better preci- 
sion of the measurements is needed in order to clarify whether the MNBD and 
EMNBD parametrizations are simply the successive approximations to the 
genuine multiplicity distribution, given by Nature or no additional iterations 
is needed. 

In conclusion, it is shown that the EMNBD describes negatively charged 
particle multiplicity distributions in e^e" annihilation into hadrons as well 
as the MNBD. The energy dependence of the parameter k, assumed to be the 
number of particle production sources, is similar for both parametrizations. 
The energy independence of the EMNBD parameters ei and £2 supports the 
toy model proposed in [1,3,4]. Better precision in multiplicity measurements 
is needed to discriminate between the MNBD and the EMNBD. 
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Figure Captions 



Fig.l The energy dependence of the parameters ei and £2 obtained from 
the EMNBD fits to the negative charged particle muhiphcity distributions 
in the e+e" annihilation. 

Fig. 2 The energy dependence of the parameter k for the EMNBD fits 
compared to the similar dependence for the MNBD fits. 



Table 1: Results of the EMNBD fits to the negatively charged particle mul- 
tiplicity distributions. The last column gives the x^ values for the MNBD 
fits taken from [4,5]. 



Experiment 


v^ (GeV) 


k 


Si 


^2 


xVndf 


X^IMNBD) 


MARKII[8] 


3.0 

4.0 

4.8 
7.4 


1±0.67 
2±0.47 
2±0.42 

q+3.8 
<J-1.4 


0.462±0.049 
0.612±0.029 
0.622±0.019 
0.544±0.074 


0.510±0.031 
0.211±0.029 
0.281d=0.018 
0.159±0.069 


5.6/2 
3.1/3 

8.2/3 
5.1/4 


3.3 

7.6 
14 
5.2 


ARGUS [9] 


9.36 

10.58 


O+0.86 
"J-O.IS 
/I+2.76 
4-0.94 


0.559±0.036 
0.596±0.119 


0.339±0.033 
0.308±0.102 


2.6/7 
1.3/7 


5.1 
1.6 


TASSO[10] 


14.0 
22.0 
34.8 
43.6 


/I +0.42 
^-0.51 
p:+0.58 
'J-0.65 
^+0.34 
O-l.Ol 
0+2.13 
°-1.05 


0.664±0.024 
0.655±0.029 
0.653±0.018 
0.567±0.030 


0.165±0.021 
0.141±0.026 
0.111±0.018 
0.086±0.034 


5.8/11 
5.2/12 
7.8/16 
12.7/17 


12 
8.2 
16 
13 


HRS[11] 


29 


fi+0.90 
O-0.93 


0.643±0.025 


0.141±0.23 


6.5/12 


7.8 


AMY [12] 


50 
52 
55 
56 

57 

60 

60.8 

61.4 


f^+1.21 
0-0.23 
7+0.69 
'-0.38 

O-0.24 
7+1.54 
'-0.92 
e+1.10 
O-0.86 
7+1.84 
'-1.20 
7+0.87 
'-0.81 
f;+0.88 
"J-0.20 


0.665±0.079 
0.724±0.143 
0.680±0.079 
0.647±0.061 
0.670±0.077 
0.653±0.092 
0.648±0.053 
0.691±0.073 


0.184±0.077 
0.041±0.181 
0.170±0.081 
0.136±0.063 
0.181±0.079 
0.120±0.101 
0.167±0.053 
0.171±0.073 


2.1/17 

6.1/17 

3.8/17 

11.4/17 

6.6/17 

6.1/18 

15.8/18 

9.1/18 


2.3 
6.1 
4.6 
12 
7.5 
6.4 
16 
10 


ALEPH[16] 


91.2 


0+1.16 
0-1.46 


0.657±0.159 


0.039±0.225 


10.9/21 


11 


ALEPH[17] 


91.2 


7+1.13 
'-2.00 


0.736±0.356 


0.035±0.444 


3.8/24 


3.8 


DELPHI[13] 


91.2 


7+0.55 
'-1.74 


0.754±0.203 


0.334±0.252 


18.6/19 


30 


DELPHI[14] 


91.2 


6±0.23 


0.639±0.022 


0.228±0.080 


94.2/23 




L3[18] 


91.2 


7+0.75 
'-0.29 


0.696±0.107 


0.074±0.119 


14.0/21 


14 


OPAL [15] 


91.2 


7+0.58 
'-0.18 


0.678±0.065 


0.101±0.066 


4.2/23 


5 


OPAL [19] 


133 


p;+1.75 
'3-0.35 


0.433±0.150 


0.493±0.113 


3.7/21 


5.3 
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